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Abstract. The physics of charm has become one of the best laboratories 
exposing the limitations of the naive constituent quark model and also giving 
hints into a more mature description of meson spectroscopy, beyond the simple 
quark-antiquark configurations. In this talk we review some recent studies of 
multiquark components in the charm sector and discuss in particular exotic 
and non-exotic four-quark systems, both with pairwise and many-body forces. 

More than thirty years after the so-called November revolution [I], heavy 
hadron spectroscopy remains a challenge. The formerly comfortable world of 
heavy mesons is shaken by new results [2]. This started in 2003 with the discovery 
of the D* (2317) and -D s i(2460) mesons in the open-charm sector. These positive- 
parity states have masses lighter than expected from quark models, and also 
smaller widths. Out of the many proposed explanations, the unquenching of the 
naive quark model has been successful |3|. When a (qq) pair occurs in a P-wave 
but can couple to hadron pairs in 5-wave, the latter configuration distorts the 
(qq) picture. Therefore, the + and 1 + (cs) states predicted above the DK(D*K) 
thresholds couple to the continuum. This mixes meson-meson components in the 
wave function, an idea advocated long ago to explain the spectrum and properties 
of light-scalar mesons [4] . 

This possibility of (csnn) (n stands for a light quark) components in D* has 
open the discussion about the presence of compact (cenn) four-quark states in 
the charmonium spectroscopy. Some states recently found in the hidden-charm 
sector may fit in the simple quark-model description as (cc) pairs (e.g., A (3940), 
y(3940), and Z(3940) as radially excited \c0, Xd, and Xci), but others appear to 
be more elusive, in particular A(3872), Z(4430) + , and y(4260). The debate on 
the nature of these states is open, with special emphasis on the A (3872). Since it 



'Article based on the presentations by J. Vijande and J.-M. Richard at the Fifth Workshop on 
Critical Stability, Erice, Sicily, Received February 2, 2009; Accepted February 9, 2009. 



2 



Four-quark stability 



was first reported by Belle in 2003 [5j, it has gradually become the flagship of the 
new armada of states whose properties make their identification as traditional 
(cc) states unlikely. An average mass of 3871.2 ± 0.5 MeV and a narrow width of 
less than 2.3 MeV have been reported for the X(3872). Note the vicinity of this 
state to the D°D*° threshold, M(D°D*°) = 3871.2 ± 1.2 MeV. With respect 
to the V(3872) quantum numbers, although some caution is still required until 
better statistic is obtained [6], an isoscalar J PC = 1 ++ state seems to be the 
best candidate to describe the properties of the X(3872). 

Another hot sector, at least for theorists, includes the (ccnn) states, which are 
manifestly exotic with charm 2 and baryon number 0. Should they lie below the 
threshold for dissociation into two ordinary hadrons, they would be narrow and 
show up clearly in the experimental spectrum. There are already estimates of the 
production rates indicating they could be produced and detected at present (and 
future) experimental facilities [7j. The stability of such (QQqq) states has been 
discussed since the early 80s [8] , and there is a consensus that stability is reached 
when the mass ratio M(Q)/m(q) becomes large enough. See, e.g., [9] for Refs. 
This effect is also found in QCD sum rules |10j . This improved binding when 
M/m increases is due to the same mechanism by which the hydrogen molecule 
(p,p, e~, e~) is much more bound than the positronium molecule (e + , e + , e~, e~). 
What matters is not the Coulomb character of the potential, but its property 
to remain identical when the masses change. In quark physics, this property 
is named flavour independence. It is reasonably well satisfied, with departures 
mainly due to spin-dependent corrections. 

The question is whether stability is already possible for (ccnn) or requires 
heavier quarks. In Ref. [9], a marginal binding was found for a specific potential 
for which earlier studies found no binding. This illustrates how difficult are such 
four-body calculations. 

In another recent investigation, the four-body Schrodinger equation has been 
solved accurately using the hyper spherical harmonic (HH) formalism |11| . with 
two standard quark models containing a linear confinement supplemented by a 
Fermi-Breit one-gluon exchange interaction (BCN), and also boson exchanges 
between the light quarks (CQC). The model parameters were tuned in the meson 
and baryon spectra. The results are given in Table [IJ indicating the quantum 
numbers of the state studied, the maximum value of the grand angular momen- 
tum used in the HH expansion, K m , and the energy difference between the mass 
of the four-quark state, E^ q , and that of the lowest two- meson threshold calcu- 
lated with the same potential model, Ae- For the {ccnn) system we have also 
calculated the radius of the four-quark state, R^q, and its ratio to the sum of the 
radii of the lowest two-meson threshold, Ar. 

Besides trying to unravel the possible existence of bound (ccnn) and (ccnn.) 
states one should aspire to understand whether it is possible to differentiate 
between compact and molecular states. A molecular state may be understood 
as a four-quark state containing a single physical two-meson component, i.e., a 
unique singlet-singlet component in the colour wave function with well-defined 
spin and isospin quantum numbers. One could expect these states not being 
deeply bound and therefore having a size of the order of the two-meson system, 
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Table 1. (ccnn) (left) and (ccnn) (right) results. 
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i.e., Ar ~ 1. Opposite to that, a compact state may be characterized by its 
involved structure on the colour space, its wave function containing different 
singlet-singlet components with non negligible probabilities. One would expect 
such states would be smaller than typical two-meson systems, i.e., Aft < 1. Let 
us notice that while Ar > 1 but finite would correspond to a meson-meson 
molecule Ar K —^° qo would represent an unbound threshold. 

As can be seen in Table [1] (left), in the case of the (ccnn) there appear no 
bound states for any set of quantum numbers, including the suggested assignment 
for the X(3872). Independently of the quark-quark interaction and the quantum 
numbers considered, the system evolves to a well separated two-meson state. 
This is clearly seen in the energy, approaching the threshold made of two free 
mesons, and also in the probabilities of the different colour components of the 
wave function and in the radius [llj . Thus, in any manner one can claim for the 
existence of a bound state for the (ccnn) system. 

A completely different behaviour is observed in Table [U (right). Here, there 
are some particular quantum numbers where the energy is quickly stabilized 
below the theoretical threshold. Of particular interest is the 1+ ccnn state, whose 
existence was predicted more than twenty years ago |12j . There is a remarkable 
agreement on the existence of an isoscalar J p = 1+ ccnn bound state using both 
BCN and CQC models, if not in its properties. For the CQC model the predicted 
binding energy is large, — 76 MeV, Ar < 1, and a very involved structure of 
its wave function (the DD* component of its wave function only accounts for 
the 50% of the total probability) what would fit into compact state. Opposite to 
that, the BCN model predicts a rather small binding, —7 MeV, and Ar is larger 
than 1, although finite. This state would naturally correspond to a meson- meson 
molecule. 

Concerning the other two states that are below threshold in Table Q] a more 
careful analysis is required. Two-meson thresholds must be determined assuming 
quantum number conservation within exactly the same scheme used in the four- 
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quark calculation. Dealing with strongly interacting particles, the two-meson 
states should have well defined total angular momentum, parity, and a prop- 
erly symmetrized wave function if two identical mesons are considered (coupled 
scheme). When noncentral forces are not taken into account, orbital angular mo- 
mentum and total spin are also good quantum numbers (uncoupled scheme) . We 
would like to emphasize that although we use central forces in our calculation 
the coupled scheme is the relevant one for observations, since a small non-central 
component in the potential is enough to produce a sizeable effect on the width 
of a state. These state are below the thresholds given by the uncoupled scheme 
but above the ones given within the coupled scheme what discard these quantum 
numbers as promising candidates for being observed experimentally. 

Binding increases for larger M/m, but in the (bbnn) sector, there is no pro- 
liferation of bound states. We have studied all ground states of (bbnn) using the 
same interacting potentials as in the double-charm case. Only four bound states 
have been found, with quantum numbers J P (I) = l + (0) , + (0), 3~(1), and 
1~(0). The first three ones correspond to compact states. 

Now, one could question the validity of the potential models used in these 
estimates, or more precisely, of the extrapolation from mesons to baryons, and 
then to multiquark states. For the short-range terms, in particular one-gluon 
exchange, the additive rule 

i<j 

is justified. Here v(r) is the quark-antiquark potential governing mesons, and X| 
is the colour generator. This is the non-Abelian version of the 1/r — > ^ QiQj/ r ij 
rule in atomic physics. 

The confining part, however, is hardly of pairwise character. Several authors 
have proposed that the linearly rising potential or of mesons (a is the string 
tension) is generalised as 

V = cr min(di + d 2 + d 3 ) , (2) 

where di is the distance from the i th quark to a junction whose location is opti- 
mised, exactly as in the famous problem of Fermat and Torricelli. Unfortunately, 
the potential (J2J) differs little from the empirical ansatz ([T]) which here reduces 
to a(ri2 +V23 +r3i)/2. Hence baryon spectroscopy cannot probe the three-body 
character of confinement. 

In the case of two quarks and two antiquarks, the confining potential reads 

V A = mm(V f , V.) , (3) 

given by the minimum of a flip-flop potential Vf and a Steiner-tree potential V s , 
sometimes named "butterfly" (see Fig. [TJ. In Vf, each gluon flux goes from a 
quark to an antiquark. The second term corresponds to a minimal Steiner tree, 
with four terminals and two Steiner points. It is remarkable that this poten- 
tial, which is supported by lattice QCD [13] is more attractive than the additive 
ansatz. This is illustrated in Ref. [Hj, where the four-body problem is solved 
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Figure 1. String model for four quarks: flip-flop (left) and Steiner-tree (right), an alternative 
configuration that is favoured when the quarks (full disks) are well separated from the antiquarks 
(open circles). 



Table 2. Four-quark variational energy E4 of QQqq for the different confinement models (V/ 
stands for the flip-flop interaction, V s for the Steiner-tree potential, and V4 = min(V/, V s )), 
compared to its threshold, and variational energy of QQqq with the flip-flop model Vf, 
compared to its threshold T4 as a function of the mass ratio. 



M/m 




E 4 




T 4 




T' 






V s 


Va 




Vf 




1 


4.644 


5.886 


4.639 


4.676 


4.644 


4.676 


2 


4.211 


5.300 


4.206 


4.248 


4.313 


4.194 


3 


4.037 


5.031 


4.032 


4.086 


4.193 


3.959 


4 


3.941 


4.868 


3.936 


3.998 


4.117 


3.811 


5 


3.880 


4.754 


3.873 


3.942 


4.060 


3.705 



with this confining term alone without short-range corrections. The results are 
displayed in Table [2 This four-body calculation is rather involved, as the po- 
tential at each point is obtained by a minimisation over several parameters. See 
Ref. [2] for technical details about the models and the numerical techniques 
used. 

The results for the configurations (QQqq) and (QQqq) are shown in Table [2] 
as function of the heavy-to-light mass ratio. Clearly, as M/m increases, a deeper 
binding is obtained for the flavour-exotic (QQqq) system. For the hidden-flavour 
(QQqq), however, the stability deteriorates, becoming unbound for M/m > 1.2. 

More recently, the stability in this model has been demonstrated rigorously 
in the limit of very large M/m. The first step is to show that 

V 4 /a < ^ (1*1 + |v|) + \z\ , (4) 

in terms of the Jacobi variables, x = — Ti, y = r± — ^3 and z = (r% + r 4 — 
r i + r 2)/2, so that the Hamiltonian describing the relative motion is bounded 

Ml ml 4/x 

(fj, is the quark-antiquark reduced mass), which is exactly solvable for its ground 
state and gives binding for large M/m. Details will be published shortly [15]. 

To conclude, let us stress again the important difference between the two 
physical systems which have been considered. While for the (ccnn), there are 
two allowed physical decay channels, (cc) + (nn) and (en) + (cn), for the (ccnn) 
only one physical system contains the possible final states, (cn) + (cn). Therefore, 
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a (ccnn) four-quark state will hardly present bound states, because the system 
will reorder itself to become the lightest two-meson state, either (cc) + (nn) or 
(cn)+(cn). In other words, if the attraction is provided by the interaction between 
particles i and j, it does also contribute to the asymptotic two- meson state. This 
does not happen for the {ccnn) if the interaction between, for example, the two 
quarks is strongly attractive. In this case there is no asymptotic two-meson state 
including such attraction, and therefore the system might bind. 

Once all possible (ccnn), (bbnh) and (ccnn) quantum numbers have been 
exhausted very few alternatives remain. If additional bound four-quark states or 
higher configuration are experimentally found, then other mechanisms should be 
at work, for instance based on diquarks [4} 116} ITT]. 
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